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Abstract 

We study the radiative decay of T into a scalar glueball T ^Gs using 
QCD factorization. We find that for this process the non-perturbative effects 
can be factorized into a matrix element well defined in non-relativistic QCD 
(NRQCD) and the gluon distribution amplitude. The same NRQCD matrix 
element appears also in leptonic decay of T and therefore can be determined 
from data. In the asymptotic limit the gluon distribution amplitude is known 
up to a normalization constant. Using a QCD sum-rule calculation for the 
normalization constant, we obtain Br[T jGs) to be in the range (1 ~ 
2) X 10^^. We also discuss some of the implications for T — > 7/j decays. 
Near future data from CLEO-III can provide crucial information about scalar 
glueball properties. 
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1 



Typeset using REVT^ 



1. Introduction 

The existence of glueballs are natural predictions of QCD. Some of the low lying states are 
0"^+, ^, 1^^ and 2++ with the lowest mass eigenstate to be 0"*""'' in the range of 1.5 ~ 1.7 
GeV from theoretical calculations |jl]]. There are indications that /o(1370), /o(1500) and 
/o(1710) contain substantial scalar glueball content. For the search of glueballs, decays of 
quarkonia are well suited processes because the decays are mediated by gluons. Among 
these decays, two-body radiative decays are ideal places to study this subject, because there 
is no complication of interactions between light hadrons. Radiative decays of T have been 
studied before, in particularly by CLEO recently. With the current data sample, there 
are already several observations of radiative decay of T into mesons. Among them only a 
few with good precisions, such as the decay T — > 7/2(1270), T — > 7/o(1710) •yKK), 
while the others have large errors M- About 4 fb"^ bb resonance data are planned to be 
taken at CLEO-III in the year prior to conversion to low energy operation (CLEO-C) [Q. 
This will produce the largest data sample of T in the world. More radiative decay modes of 
T may be observed. Combining experimental data in the near future and theoretical results 
glueball properties can be studied in details. 

In this paper we carry out a theoretical study of the radiative decay of T into a scalar 
glueball by using QCD factorization. We find that the non-perturbative effects can be 
factorized into a matrix element well defined in non-relativistic QCD (NRQCD), and the 
gluon distribution amplitude. The NRQCD matrix element can be determined from lep- 
tonic T decays. The asymptotic form of the gluon distribution amplitude is known in 
QCD up to a normalization constant. Using a QCD sum rule calculation for this con- 
stant, the branching ratio Br{T — 'yGs) is predicted to be in the range of (1 ~ 2) x 10^^. 
Combining this result with experimental data, we find that none of the candidate scalar 
glueballs /o(1500) and /o(1710) can be a pure glueball. Existing information on glue- 
ball mixing allow us to predict the branching ratios for several radiative decays, such as 
T — > 7/0(1370,1500,1710) •yKK^mr). A mixing pattern suggested in the literature is 
shown to be in confiict with data. Near future experimental data from CLEO-III will pro- 
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vide crucial information about scalar glueball properties. 



2. QCD factorization of T ^ 

It is known that properties of T can be well described with non-relativistic 
QCD(NRQCD) [Q. The decay of T ^ 76*5 can be thought of as a free bb quark pair 
first freed from T with a probability which is characterized by matrix elements defined in 
NRQCD, this pair of quarks decays into a photon and gluons, and then the gluons subse- 
quently converted into a scalar glueball. In the heavy quark limit 00, the glueball 
has a large momentum, this allows a twist-expansion to describe the conversion. Also, the 
gluons are hard and perturbative QCD can be applied for the decay of the bb pair into a 
photon and gluons. This implies that the decay width can be factorized. In the real world, 
the b-quark mass is 5 GeV and a scalar glueball has a mass around 1.5 GeV as suggested 
by lattice QCD simulations 0. This may lead to a question if the twist expansion is ap- 
plicable. For the radiative decay of T, the glueball has a momentum of order of rrib- The 
twist expansion means a collinear expansion of momenta of gluons in the glueball, compo- 
nents of these momenta have the order of {0{k'^), 0{{k~), 0{Aqcd), 0{Aqcd)), where k is 
the momentum of the glueball. Here we used the light-cone coordinate system. Hence the 
expansion parameters are 



k tuq a 
F=M?^°-°2' -k 



0.02, -f^~0.1, (1) 



where is the mass of the glueball and we have taken Aqcd ~ 500MeV. In the above 
estimation we have used the fact that the probability for the conversion of gluons into a 
glueball is suppressed if the "+" component of the momentum of a gluon is very small. We 
see that the relevant expansion parameters are small, therefore twist expansion is expected 
to be a good one. We note that the same approximation may not be applied to J/ ip system 
because in this case the relevant expansion parameters are not small. We now provide some 
details of the calculations. 

The leading Feynmann diagrams for T 'yGg are from bb annihilation into two gluons 



and a photon. The basic formahsm for such calculations have been developed in Ref. |0 and 
has been used in the case of T ^ iviv') to obtained consistent result with experimental 
data 1^. With appropriate modifications we can obtain the S-matrix for T '-/Gg decay. 
It is given by 

{jGs\S\T)= -i^eQ.gX J d''xd^yd''zd''x,d%e'^-^Gs\G';^{x)Gt{y)\0) 

mixMy^)\T) ■ Ml^'''''\x,y,x^,y,,z), (2) 

where Mj^J^^'"'^ is a known function from evaluation of the Feynmann diagrams, i and j stand 
for Dirac- and color indices, a and b is the color indices of gluon. e* is the polarization vector 
of the photon and Qb = —1/3 is the b quark electric charge. Since b quark is heavy and 
moves with small velocity v, one can expand the Dirac fields in NRQCD fields: 

{0\b,{x)kiy)\T) = -i(P_,Vp_),^.(O|xV0|T)e-^-(-+^) + Oiv'), (3) 

where xH"^) is the NRQCD field for b{b) quark and P± = (1 ± 7*^)/2. The b is almost at 
rest, then = {nib, 0, 0, 0) with rUb being the b quark pole mass. 
From the above we obtain the decay amplitude for T '-/Gg as 

T='-^{0\xh-aij\T) fdz—^Ts{z), (4) 
6 JO z[l — z) 

the decay width then reads: 
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(5) 

In the above Ts is the gluon distribution amplitude of Gg and is given by 

Here we have used a gauge with G~^ = such that the gauge link between the field strength 
operators vanish. This distribution characterizes basically how two gluons are converted 
into Gg, where one of the two gluons has the momentum {zk~^ , 0, Ot)- 
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9m^ 



n'Qlaal{T\0,CS,)\T) 



dz- 



z(l-z 



-Ts{z) 



In the above equations, the matrix element (T|0i(^5'i)|T) is defined in NRQCD contains 
the bound state effect of b-quarks in T and can be extracted from leptonic T —>■ 
decay. A prediction can be made for T '-/Gs if the distribution amphtude is known. 

The distribution amphtude can be written as 

J^,{z) = fj{z), with fdzf{z) = l. (7) 

Jo 

where f{z) is a dimensionless function and its asymptotic form is: 

f{z)=?,Qz\l-zf. (8) 

With the asymptotic form in Eq.(6) we have: 

^ r(T^7 + a) ^ 257ra| 

r(T ^ (+(r) 3« 'ml- ^ ' 

In the above we have used the fact that both T jCg and T — > l^l^ are proportional to 
(T|Oi(35i)|T). 

The use of the asymptotic form for f{z) may introduce some errors, because the scale 
fi here is actually rUb, not fi oo. However, with a large ruh one may expect that it can 
provide a good order of magnitude estimate with the asymptotic form. We will use Eq. @ 
later for our numerical discussions. 

We note that at this stage the state Gg can be any particle with the same quantum num- 
ber as G°''~^^G1^'^ , i.e., J^'" = 0^'^. The normalization constant fg depends on the properties 
of the specific particle. In order to obtain the branching ratio of Gg as a scalar glueball, we 
have to evaluate fs with Gg specified to be the scalar glueball. In the following we provide 
an estimate based on QCD sum rule calculation. 

3. QCD sum rule calculation of the normalization constant 

The constant fg has dimension one in mass and is related to the the product of local 
operator: 

{Ggik)\G^PG%\0) = Umlg^- + /.PF. (10) 



The fact that the same fs appears in Eq.(9) and Eq.(12) can be checked by integrating over 
z on the both sides of Eq.(6). 

The basic idea of the QCD sum rule calculation for our estimate is to consider the two 
point correlator 

n^.,MV'(Q') = / d^xe'''-^i{0\TG,o.G:{x),G,,pG^,,mo) 

+Sl,'u'^s2{Q') + Sl^,^,Uss{Q'), (11) 

for a region of Q in which one can incorporate the asymptotic freedom property of QCD via 
the operator product expansion (OPE), and then relate it to the hadronic matrix elements 
via the dispersion relation. The tensors in Eq. (|n|) are defined as 

t t t t t t 

Vf^ut^'u' = 9lf,'(iuqu' + 9L'%%' + 9lu'(iuqi,' + 9if,'(itiq,y' 

^liufi'u' = 9liu9l.'u'j ^l.ufi'u' = 9liuQiJ''Qv' + 9li'u'Qfj.Qv, S^ufi'u' = Qfj-QuQfi'Qu' , (12) 

where gj^^ = g^^ - q^Qy/q^. The corresponding terms TIt{Q'^), TiyiQ'^), ^si{Q^), ^S2{Q'^) 
and TisiiQ"^) are from the contributions of 2"''+, 1 ^ and 0"^"^ states respectively. 
In a deep Euclidean region = —q^ » Aqcd, they can be expanded as 

n.(g') = C^{Q')I + Cl{Q')as{G,yGn + Cf{Q'){gsr"^G^ ^^G" ^^G^ + ■ ■ ■ , (13) 

where Cf are Wilson coefficients which need to be determined later. 

On the other hand, the correlator in Eq.(^) can be saturated by all possible resonances 
and continuum. We have 

Jmn^,,^v'(Q') = 5](0|G^„Ge|i?)(i?|G^';3G(^.|0)7r5(g2 + rn^) + continuum, (14) 

R 

where the sum on R are for all possible resonances. The term {\G ^aG'^\R) {R\G i3G^,\Q) 
in the above equation contains the information of /o and fg when R is the scalar glueball. 
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The T and V tensors are not related to fs- They are irrelevant to our calculations. The 
functions n5'i ,52,S3 contain linear combinations of /o and fg. QCD sum rule calculations 
for {Gs\G^'^G^iy\0) = (4/0 + fs)TTT'G has been carried out before 0. Therefore if one of the 
nsi,s2,S3 is know, one can obtain fg. From eq. (10) and the tensor structure of eq. (11), 
we find that 1153 is proportional to (/o + fs)"^- Therefore the study of Uss is sufficient for 
our purpose of determining fs- Ilsi,2 also contain information about /o and fs- However the 
non-perturbative contributions for them begin at the level of dimension-8 operators. The 
results obtained are not as reliable as the one from Uss which has a lower dimension. We 
now concentrate on Uss- 

There may be several bound states with the same quantum numbers to include in the 
QCD sum rule calculation, such as a pure scalar glueball, quark bound states and higher 
excited states. The contributions from higher excited states are suppressed upon the use of 
Borel transformation which is discussed in the below. For the quark bound states, OZI rule 
implies that the conversion of bound quark state into a scalar glueball is suppressed compared 
with the conversion of two gluon into a scalar glueball |T^, perturbatively suppressed by 



powers in a^. If this is indeed true, the corresponding fg parameters for quark bound 
states will be smaller than pure glueball state. We will work with this approximation in the 
following discussions. To be consistent with our previous expansion, we again work to order 



ttg- To this order, using the method in Ref. |11|, we find 



^ssiQ') = ^1^^ + ^(^^M.^"'^) + ^(r'^G'^ -,G^ (15) 



The correlator in Eq. ([T5| ) obtained by using OPE is related to Eq. (|T4D via the standard 
dispersion relation 



n,.,,v(Q^) = - / ^^'^r: ^ (16) 



n Jo s + Q'^ 

In practice one may only include ground states in the calculation. In order to reduce uncer- 
tainty due to higher excited states and also continuum states, we apply Borel transformation 
and obtain 



BIlssiQ') = fj dse-^/'^'pssis), (17) 
where ^53(5) = (l/7r)/mn53(— s), and 

BUiQ') = lim -^(Q'r(--±-)-UiQ'), (18) 



Here one also needs to have the hmit Q'^/n = = constant. 

In our numerical calculation we have varied sq in the range of 3 — 6GeV^, and found that 
the uncertainty is around 10 percent. The parameters determined are reasonably stable. 

We obtain the range for fg as 

= (100 ~ 130)MeV, (19) 

with /o = 190 MeV and /, = 100 MeV for mo++ = 1.5 GeV, and /o = 130 MeV and 
fs = 130 for mo++ = 1.7 GeV. In obtaining the above result, we have re-evaluated /o also 
using the same parameters. The input parameters used are [1^ «s(yu) = 4:7r/9\n{fi'^ / Aq^j^), 
Aqcd = 0.25GeV, /i = M, {asG^^G"^) = 0.06 ± 0.02GeV\ and c/,(r*"G" '^^G^ "^G^ = 
0.27GeV^{asG''^G^^). 

For consistency, we also calculated the glueball masses. We find that for the 0++ state 
the mass is 1.5 ~ 1.7 GeV, and for 2++ the mass is 2.0 ~ 2.2 GeV. These are in agreement 
with other calculations f^. 

If the scalar glueball is a pure one, using the above results we obtain the branching ratio 
for T — i> '-fGs to be in the range 

Br{T -fGs) = (1 ~ 2) X 10-^ (20) 

with a larger branching ratio for a larger glueball mass up to 1.7 GeV. Here we have used 
as = 0.18 which is the typical value for in the energy range of the decay. We obtain 
a large branching ratio for T 'yGg- We would like to point out that considering several 
uncertainties, the assumptions of factorization and single pure glueball state in the QCD 
sum rule calculation, the above numbers should be used as an order of magnitude estimate. 
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4. Discussions of phenomenological implications 

Experimental measurement of T ^ '-yGs may be non-trivial. One has to rely on the 
decay products of glueballs. There are several ways the glueball can decay with reasonably 
large branching ratios, Gs —>■ KK or Gs to multi-pions. As mentioned earlier that there 
are several candidates for scalar glueball, the /(1370), /o(1500) and /o(1710). Decays of 
T 7/0 (^) 'y{KK or multi-pions) can provide important information. 

Experimentally there is only an upper bound of Br(T 7/o(1710) — > jKR) < 
2.6 X 10^'' at 90% C.L.. If /o(1710) is a pure glueball, experimental measurement [H of 
5r(/o(1710) ^ KK) = O.SStHg 1 would imply Br{T 7/o(1710) ^ -fKK) to be in 
the range of (0.4 ~ 1.0) x 10^^ which seems to indicate that /o(1710) may not be a pure 
glueball. At present it can not rule out the possibility that one of the /o(1370) or /o(1500) 
being a pure glueball state. Data also allow certain mixing among glueball state and other 
quark bound states. 

Theoretical calculation of the mixings among glueball and quark bound states is a very 
difficult task. There is no reliable theoretical calculation. Lattice calculations may eventually 
give accurate predictions for the mixing parameters. At present there are some phenomeno- 
logical studies of glueball mixings. We now study some implications of the branching ratio 
for the radiative decay of a T into a pure scalar glueball obtained in the previous section on 
a mixing pattern suggested in Ref. ||T3[| . 

An analysis combining other experimental data in Ref. showed that the three O^"*" 



states /o(1370), /o(1500) and /o(1710) all contain substantial glueball content. Ref. |T3 



obtained a mixing matrix of physical states in terms of pure glueball and other quark bound 



states to be 13 



Jil Ji2 JSi 

/o(1710) 0.39 ±0.03 0.91 ±0.02 0.15 ±0.02 
^ ' (21) 

/o(1500) -0.65 ±0.04 0.33 ± 0.04 -0.70 ± 0.07 

/o(1370) -0.69 ±0.07 0.15 ±0.01 0.70 ± 0.07 
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where the states Gg, S = \ss) and = \uu + dd)/\/2 are the pure glueball and quark bound 
states, f^'' indicate the amphtude of glueball Gg in the three physical fo{i) states. 

Because of the mixing, when applying our calculations to radiative decay of T into a 
physical state which is not a purely gluonic state the parameters will be modified. If the 
mixing parameter is known one can obtain the Rg ratios for T — 7/o(1370, 1500, 1710) as 

i?.(T-7/(^)) = ^-^l/-f^P- (22) 

T ifii) niay also result from T decays into a 7 and 5, quark bound states. However, 
these processes are suppressed by a^. 

Using the mixing amplitudes in Eq.(^), one obtains the branching ratios of, 
T ^ 7/0(1370, 1500, 1710) to be in the ranges, (4.8 ~ 9.6,4.2 ~ 8.4,1.5 - 
3.0) X 10-^ Combining the branching ratios of /o(1370, 1500, 1710) ^ KK{'mi)) = 
(0.38lH^(0.039l°:°°I),0.044l°:°il(0.454l°:l°^),0.35l°:li(0.26l°^^ i, we obtain: 

Br{T 7/0(1710) ^ -fKK) ^ 0.6 ~ 1.2; Br{T 7/o(1710) ^ jTin) ^ 0.06 ~ 0.12; 
Br{T ^ 7/0(1500) -fKK) ^ 0.2 ~ 0.4; Br{T 7/o(1500) ^ -fnn) ^ 1.9 ~ 3.8; 
Br{T 7/0(1370) ^ -fKK) ^ 1.7 ~ 3.4; Br{T 7/o(1370) -> -fnn) ^ 1.2 ~ 2.4. 

In the above the branching ratios are in unit 10~^. The branching ratios predicted above 
can provide further test for QCD factorization. Future experimental data from CLEO III 
will provide us with important information. 

To summarize, we have estimated the branching ratio of T —>■'-/ + Gs with Gs as a 
glueball. Our result shows that /o(1710) may not be consistent with the assumption that 
it is a pure glueball, but can not rule out the possibility that one of the /(1370), and 
/o(1500) being a pure glueball state. We also predicted several T — > ■yKK^nn) branching 
ratios using a phenomenological glueball mixing pattern which can provide further tests for 
QCD factorization calculations and glueball mixing. To have a better understanding of the 
situation, we have to rely on future improved experimental data. Fortunately CLEO-III will 
provide us with more data in the near future. We have a good chance to understand the 
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properties of scalar glueball. We strongly encourage our experimental colleagues to carry 
out the study of radiative decay of T into a scalar glueball. 
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